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SUMMARY 

The dynamic and composite nature of propagation impairments that are incurred on 
earth-space communications links at frequencies in and above the 30/20 GHz Ka band, 
i.e., rain attenuation, cloud and/or clear air scintillation, etc., combined with the need to 
counter such degradations after the small link margins (typical of such frequency bands) 
have been exceeded, necessitate the use of dynamic statistical identification and prediction 
processing of the fading signal in order to optimally estimate and predict the levels of each 
of the deleterious attenuation components. Such requirements are being met in NASA's 
Advanced Communications Technology Satellite (ACTS) Project by the implementation of 
optimal processing schemes derived through the use of the ACTS Rain Attenuation 
Prediction Model and nonlinear Markov filtering theory. The ACTS Rain Attenuation 
Prediction Model discerns climatological variations on the order of 0.5° in latitude and 
longitude in the continental U.S. The static portion of the model gives precise availability 
predictions for the "spot beam" links of ACTS. However, the structure of the dynamic 
portion of the model, which yields performance parameters such as fade duration 
probabilities, etc., is isomorphic to the state-variable approach of stochastic control theory 
and, as shown in this paper, is amenable to the design of such statistical fade processing 
schemes which can be made specific to the particular climatological location at which they 
are employed. 



I. INTRODUCTION 


The development of the ACTS Rain Attenuation Prediction Model [1,2,3], another 
in a variety of such models that now exists, was in response to several system design and 
performance requirements of satellite systems in and above 30/20 GHz (the Ka band); Such 
a satellite system is NASA's Advanced Communication Technology Satellite (ACTS) 
Project after which the model is named. Two major requirements are the ability to 1) 
estimate the level of a communication signal fade due to rain on a satellite link using, in the 
general case, attenuation measurements that may be corrupted with random as well as 
systematic measurement errors (defined here as identification ) and 2) predict what such 
attenuation levels will prevail a short time into the future so as to forewarn the need for the 
deployment of a fade countermeasure (defined here as prediction ). For example, a situ- 
ation may exist for the user of a small remote satellite terminal to employ the received com- 
munications channel, with all its attendant power fluctuations due to modulation, etc., as 
the source of measurement of link attenuation due to rain so as to drive some fade mitiga- 
tion technique, the use of which is needed beyond some pre-established fade threshold. 

Or, as in the case of ACTS, one may be operating at a frequency that is not only impaired 
by rain but also by the phenomena of clear-air and/or cloud scintillation; here, one receives 
the total fading signal and, if one is to have a reliable satellite communications link, must 
separate out the component rain and scintillation effects since each must be dealt with in a 
different manner (e.g., rain fade by power control and scintillation, if it proves to be a 
problem, by time diversity transmission). 

It is the purpose of this paper to indicate how results of the dynamic portion of the 
ACTS Rain Prediction Model can be used in the implementation of processing schemes that 
are robust to such scenarios indicated above. In particular, the simplest case is considered 
whereby one has available attenuation measurements with associated measurement errors at 
discrete points in time. From these discrete, "noisy" observations of link attenuation, it is 
required to obtain optimal estimates of not only the satellite link attenuation that corre- 
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sponded to the measurement, but also an optimal prediction of what the attenuation value 
will be at the next (future) sampling time. This problem is formed within the context of 
non-linear Markov filtering. The optimality criterion used will be the minimization of the 
least square error that exists between the predicted value of attenuation and that estimated 
from the noisy measurement process. The only fading mechanism that will be considered 
here is that due to rain; scintillation and other simultaneous signal power impairments can 
also be considered but at the expense of a much more complicated and involved exposition 
Such cases will be differed to future publications. 

II. A BRIEF OVERVIEW OF THE DYNAMICS PORTION OF THE 
ACTS RAIN ATTENUATION PREDICTION MODEL. 


From the development of the ACTS Rain Attenuation Prediction Model (in particu- 
lar, Section 5 of [3]), one models the temporal evolution of the link attenuation A(t) 
through the parameter x A (t), i.e., 


^(0 = 


In A(t ) - In A m 

G lnA 


( 1 ) 


where A m is the median of the link attenuation and <J lnA is the standard deviation of the 

logarithm of attenuation; these two parameters are specific to the particular location, fre- 
quency of operation, and geometry of the satellite link. It is the x A (t ) parameter that is 

given by, in the most general case, a multi-component Markov random process which is 
determined by a system of first order stochastic differential equations. It is, however, ex- 
pedient to reduce to a one component model with careful consideration given to the incor- 
poration of the temporal smoothing process induced by the extended propagation path. In 
particular, as was done in [3], one can consider (approximately) the one component pro- 
cess given by the single first-order differential equation 
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( 2 ) 


< ^jf = -Vs x A+J 2)^S(0 
where the random function Q t) is governed by Gaussian statistics with a zero mean (i.e., 
"white noise"): 

<£(')>= 0, <|(t)|a+T))=5(r) 

and where y s is a "smoothed" temporal parameter that is given by the solution of the tran- 
scendental equation 




+ exp 

si 

exp 


— 


1 ?sj 


rs) 


-exp (— l) 


(3) 


From geometrical considerations concerning the propagation path, characteristic rain cell 
size, and long term spatial isotropy of rain cell movement, one has the following relation 
that exists between the coefficients y^ and y 2 that appear in Eq.(3): 


>- 1/2 


72= 7i 


1 + 


2Lcos8\ 
nR r 


(4) 


where L is the total propagation path length within the potential rain region, 8 is the link 
elevation angle, and R c = 4 Km is the characteristic rain cell size as required by funda- 


mental considerations and model constraints [3]. By the same considerations, one has for 
the coefficient y 1 


2v -i 

7 j = = 0.1336 min (5) 

nR c 

where v = 14 m/s is the characteristic speed of the rain cell that is assumed in the model. 
For a typical communications link with an elevation angle ~ 40°, y s ~ 0.06 min' 1 = 0.001 

sec' 1 . 
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In what follows, the one component Markov process x A (t) given by Eq.(2), with 
contact made to the link attenuation Aft) through Eq.(l), will be used to derive an optimal 
fade detection and prediction algorithm for use with measured (real time) fade data. 

III. DEVELOPMENT OF AN OPTIMAL FADE IDENTIFICATION AND 
PREDICTION ALGORITHM 

A. The Observation and Sampling Model 

Here, the problems of fade identification and prediction as defined in the introduc- 
tion will be solved. At the outset, it is useful to consider the measurement process via the 
following observation model : 


A ob Jit) = A(t) + n (t) (6) 

where Aft ) is the actual attenuation that exists on the communications link at time t, n(t ) is 
the associated measurement uncertainty or "noise", and A a bdft ) is the observed link attenu- 
ation. The measurement noise, which can result from inaccuracies in the hardware of the 
measurement process or from an assumed measurement function such as that used in 
"frequency scaling" whereby a signal fade at one frequency is derived from a fade mea- 
surement at another frequency, is a random function of time and can also be statistically 
characterized by a white noise process, viz, 

<n(r)> = 0 

2 (7) 

(nftjn (t^y=c n Sft l - t^) 

where o n is the standard deviation of the measurement noise. Remembering Eq.(l), one 
has from Eq.(6), 

A obl l ) = Mx A ,f) +n(t ) , A(x A ,t) = A m exp ( a inA x A(0) (8) 
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for the observation model where the random process x A is governed by the stochastic dif- 
ferential equation Eq.(2). The following problem can now be defined: given the measured 
fade A 0 b s (t ) and given Eqs.(8), (7), (2), and (1), it is desired to obtain an estimate of the 
quantity x A (t ) at time t (i.e., identification) as well as the extrapolated estimate of the 
quantity x A (t + x) that should prevail at a future time t + x (i.e., prediction). Using 
Eq.(l), one can then easily obtain the corresponding values A(t) and Aft + x) . 

This is formally a problem in the optimal (with respect to a given criterion) estima- 
tion (or filtering) and prediction of a non-linear Markov continuous random process sam- 
pled in time. 

Since the random process in question is the quantity x A and it is desired to obtain 
information concerning the instantaneous value of this quantity from the measurement of 
A 0 bs, one must invariably deal with a parameter known as the a posteriori probability den- 
sity p(x A (t) I A 0 bs(t) ) governing the random process at a time t conditioned on the obser- 
vation of the quantity A 0 bs(t) . Once this probability density has been secured, one can 
easily obtain descriptions of the random process^ in the form of statistical moments such 
as the average x A (t) ( referred to as the optimal estimate minimized with respect to the 

4c 

mean-square error), the associated standard deviation ) of the optimal estimate, etc. 

Although both x(t) and A 0 b s (t) are, in general, continuous in the time variable, the 
latter is actually sampled at discrete times t i with a corresponding time interval T separating 
the samples, i.e., T = t { - fa . Thus, due to this sampling process overlaid on the contin- 
uous process A 0 b s (t) , one must also now consider additional random processes which re- 
sult from the sampling, viz, 


A ob Jit^=A ob Jit)S(t-t i ) and x A ftf) =x A (t)Sft-t i ). (9) 

One must thus amend the a posteriori probability density p(x A (t) I A 0 bs(t) ) to reflect this 
circumstance and consider p (x A (t ) , x A ft) I A obs ft )), viz, the probability density gov- 
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eming the continuous random process x A (t) and the sampled process x A (t ) conditioned 
on the value of the discrete measurement of A 0 b^t) . 

B . The Statistical Identification of The Prevailing Link Attenuation 
Level 


Now the random process x A (t) , described by Eq.(2), has associated with it a tran- 
sition (and conditional) probability density p [x A (t) | x A (/*)) giving the statistics connected 
with the evolution of the sampled value x A (tj) at a time t i to the value of its continuous 
counterpart x A (t ) a later time t > t-. Following the results in [2,3] (in particular, Eqs.(53) 
and (54) of [3]), this transition probability density is given by the Kolmogorov equation 
associated with the first order differential equation Eq.(2), i.e.. 


dp (x A (f)l x A (/,•)) 


D x ^p(x A {t )l* A (fi))], 


where the differential operator D x ^ ]f (x A )] for the Kolmogorov equation is given by 


o x Jf m] ■ M + 7 S 

OX A dx A 


In an analogous fashion, the probability density p (x A (t ),x A (t ) I A obs {t )) is a transition 
probability density governing the random processes x A (t) and x A (t) conditioned on the set 
of measured observations A obs (t ), each of which is connected with the random process in 

question via the observation model of Eq.(6) and the sampling model of Eq.(9). In this 
case, the Kolmogorov Equation of Eq.(10) must be augmented with this (non-linear) ob- 
servation model. In particular, amending the observation model of Eq.(6) with the sam- 
pling model of Eq(9), one has that the quantity 
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n(t) = A,aO-Ad A ,t) (12) 

2 

is a zero mean Gaussian random variable with variance o„ (by Eq.(7». As shown in the 

Appendix, using this fact in an extrapolated a posteriori probabilistic analysis for the evo- 
lution of x A (t ), one obtains an integrodifferential equation for P ix A it) ^ x A it) \ A obs it )) 

a modified version of the Stratonovich Equation [4], well known in Markov filtering the- 
ory, viz. 


|^/t if ) » x A (t ) I A obs ( f,' ) j 

dt 


D x jp[x A it),x A it)\A obs (r, ))] + ( x a it) ,t i ) - 

~f f 0 ( * / (O , )p[x A ' (t),x A ' it)\A obs (ti )) d x A dx A 

x p jxyj (t ) , x A it ) I A obs it j ) 


(13) 


where 


1 / \2 

<p ( x A it ) , ti ) = -\A ob £ti ) - A(x ,ti )) (1 4) 

It is important to note that Eq.(13) holds for times t within the interval ti<t<t i+] . 

Obtaining a solution from Eq.(13) commences with writing the a posteriori proba- 
bility density in the form 

p\x A it),x A it )\A obs it i ))=p(x A it)\x A itj)p{x A it)\A obs it i )) (15) 


Substituting Eq.(15) into Eq.(13), integrating the result with respect to x, and using the 
facts that 

p(x A it)\x A it ))dx A = l and 



I D x \p (*A it)\x A it ))] dx A = 0 

J -M 
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(such integrals indicate integration over all the possible values realized by x A at a fixed time 

t and where the latter relation is obtained using Eq.(l 1) and the boundary conditions that 
p (*a (f ) I x a (* )) and its first derivative vanish at x = ± °°) yields the relation 


- f 0 ( X A (f ) , ti ) P ( x A ' (t) 1 A obs (ti )) d x A 


p(x A (t)\A obs (ti)) (16) 


for tj<t< t i+} . By substituting Eq.(16) back into Eq.(13), one then re-derives Eq.(10) 
prescribing the probability density p [x A (r ) I (f)) thus making the proposed solution of 

Eq.(15) self-consistent. 

An approximate analytic expression can be obtained for Eq.(16). One notes that the 
second term within the brackets is simply the average <^<p( t £) over all possible values of 
x A . Since it is only a function of t, and if the time interval t -ti is small as compared to 
the characteristic time if variation for this average, one can neglect it and obtain a simpler 
form which has the solution 


p(x A (t)\ A obs (ti)) = N exp J <p(x A (t ' ),ti ) dt ' 


where the normalization constant N is given by 


p(x A (t i )\A ob £t i j) ( 17 ) 



f 

<p(x A (t'),ti)dt' 


p(%A(ti)\A 0 blti)) 


(18) 


Equations (17) and (18) define the "partition" representation of the a posteriori probability 
density function [5]. 
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One now admits the well established Gaussian approximation, which demands that 
the probability density p ( x A (/,• X'W' ,•)) take the form 


/ 

P $ A (U ) I A ob pi )) = A / exp 

V lx K{t i) 

\ 


2 

2K(t t ) 


(19) 


where x A (t, ) is the extrapolated estimate of the random process x A (f, ) available at time t i 
and K(t i ) is the standard deviation of this estimate about the actual value; the extrapolated 
estimate, which is obtained from the optimal estimate of the previous time interval, is em- 
ployed here since it is the only a posteriori information that one has at time r,- about the 
random process x A ( t ) to combine with the observed quantity A obs (ti ) to yield the statis- 
tics that govern the possible values for x A (t t ). For this reason, K(t) is sometimes re- 
ferred to as the standard deviation ( or , in general, the covariance) associated with the er- 
rors of the extrapolated estimate. To make the problem defined by the substitution of 
Eq.(19) into Eqs.(17) and (18) analytically amenable, one expands the function 

${*^(*1 )>*i ) that appears within the exponential functions in Eqs.(17)and (18) into a Taylor 

, 

series about the optimal estimate x A (t ) obtained at the i th sampling time t i , retaining terms 
up to the second to match that of Eq.(19). Using Eq.(9) in Eq.(14) and substituting this 
composite result into the expansion and using Eq.(19), one obtains from Eqs.(17) and (18) 


P 



exp 


a (U )(*a “ *) ~ 


I A 


obs 


«i>)= 




Ga-x) 

2<U',) 


I 


V 2no (ti ) exp 


A 

i 2 / 


( 20 ) 
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where the quantity c (f f ) is the error covariance of the optimal estimate (the latter will be 
defined shortly) given by 

&)-&(*/)) (21) 

and 


a (0 s 


d<PJS, 0 


_ 


dx 


, b(t}m 


1 d 


dx 1 


( 22 ) 


are coefficients that contain the observed sampled values A obs ( t L ); here, 
Qsi x a ( f t)> *i ) = $( x a ( 0 > t - f,) . Using Eq.(14), one then obtains 




a(ti) = 


'in A 


A (x ,f f )(a ob lti )-A(x a ,ti )) 




/ 2 \ 

C lnA 


(23) 


A(x,t i iA ob £! i )-2A$ A ,t i )) 


Using the a posteriori probability density of Eq.(20), one can now obtain an ex- 
pression for the optimal estimate x *(t-) of the value of x A (r t ) connected with the observa- 
tional data A obs (r, ). Of the many optimality criteria that can be used to define a particular 


estimate, the one selected here is the simplest, viz, the criterion that minimizes the mean 
square error which yields the relation 



X A (*i>p[xA(fi) iA obs^) d ^ 


Substituting Eq.(20) into this and performing the integration yields the following recurrsion 
relationship for the optimal estimate: 


(24) 


x/Kfd^xM + Oxpbaitd 


As mentioned earlier, the extrapolated estimate x A (*») is that obtained from the op- 
timal estimate of the previous sampling time through the use of the transition probabil- 
ity density defined by Eq.(10). The solution of this equation is well known [2,3] and is 
given by 


p(*a (r i )lx i 4 (r,_i))= 



% A (fi) * A (fi — l) 


2 1 - 0 . 


(25) 


where 0; = exp ( -y s (f ; - t,_i )). The extrapolated estimate x A (tj) at time f,- given the a 
~ * 

posteriori value x A (f,_j) = x A (t;_ j) is simply 


x A (ti> 



X A P^AiO^A^i-Oj^A 


Upon substituting Eq.(25) into this relationship and integrating, one obtains 


/ 


x A (ti) = ®iX A (ti-i) 


(26) 


which, when used in Eq.(24), demonstrates the recurrsive nature of the estimate. 

It now remains to obtain a similar expression for the standard deviation of the ex- 
trapolated estimate K(t ) . Equation (19) essentially provides a working definition for this 
quantity. In particular, one has 
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K(ti) = 

In order to relate a value of this parameter at ti to the a posteriori values obtained at time 
j, one can use Eq.(15) to write the relationship 



K(ti) = a (x A (h ) - *A (ti )) P @A 7 (ti),*A (ti- 1 ) 1 A ob^i- 1 ))dxdx' (27) 


Using the previously derived expressions of Eqs.(19) and (25) in Eq.(27) yields the re- 
currsive relationship 


K«i)= l + ( 28 ) 

* 

that relies on former values of O x (ti ) obtained from Eq.(21). 

Equations (24), (21), (26), and (28), in addition to the auxiliary relations of 
Eq.(23) collectively compose the dynamic identification part of the problem defined at the 
outset. The approximation used that allows one to write Eq.(17) as well as considerations 
of the errors incurred in the optimal estimate place upper limits on the time interval [tj. j , ti) 
between consecutive samples. However, the typical clock intervals ~ 1 second that are 
encountered in practice are easily within this limit. It now remains to obtain a prescription 
for the dynamic prediction of link fade levels into the future based on the data afforded 
from this fade identification process. 

C. The Prediction of Link Attenuation for Short Times into the Future 

From some results of the foregoing, in particular, the solution given by Eq.(25) for 
the transition probability density of the attenuation process, one can derive a relation yield- 
ing the average value of the link attenuation at a time T pred into the future. In particular, the 

problem can be succinctly defined as follows: Given the fact that an optimal estimate x A (t t ) 


1 3 



is obtained for the attenuation process at the time t p and on the hypothesis that the attenua- 
tion will increase , it is of interest to obtain the extrapolated estimate Ax A ( t t +T pred ) of the 
change in the attenuation that will occur at a later time f; + T pred . One then has for the total 
value of this combined extrapolated optimal estimate of the attenuation process 

X A (*i +T pred) - X A (0 + ( f i +T pred) (29) 


By this formulation of the problem, one has 


^A^i +T pred) = 


I 


('< + T pred) P \ X A (U +T pred) > x A (0 + A (*i +T pred) 1 X A (t) U (Ax a ) 


(30) 


where the conditional probability density is given by 


P 1 * A (U p r ed) ^ % A (*i) A ^ pred) ^ X A I 


I. 


p(x A '(t)\x A (t i )dx A / (31) 


XAitJ + AxJti+TpJ) 


To this end, substituting Eq.(25) into Eq.(31) and this result into Eq.(30) yields, 
upon a change of variables, 


A (h +Tpr e d) ~ 


f. 


erfc\ 




V2D W^d) 


\ydy - 
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where & pred = exp f-y s T pred ) and D (T pred ) =1 - <P pred . In the approximation in 
which one can assume that x A ( 1 —® pred ) ~ 0 relative to the range of integration and the 


contributions of the integrands, the integrals in Eq.(32) can be analytically evaluated and, in 
conjunction with Eq.(29), gives for the extrapolated optimal estimate of attenuation for a 
future time t, +T pred based on the estimate at time t t 


*>, = V pnd (,,) + (33) 

where 

(34> 


A limitation is placed upon the maximum value of T pre d by the approximation used 
to evaluate the integrals of Eq.(32); in particular, one has 


Tpred « 


' InA 


2Y s [lnA{t)-lnA m } 


showing that, as can be expected, the value is a function of the prevailing link attenuation 
A(t). A more rigorous analysis for this maximum value would involve considering the 
prediction error associated with Ax A (r,- +T pred ). However, since in the applications, one 
usually only needs to predict no more than two to four time intervals ahead (i.e., typically ~ 
2 to 4 seconds), this condition is easily met. 
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IV. ALGORITHMIC IMPLEMENTATION OF THE FOREGOING 


The synthesis of the relevant relations derived above is shown in Fig. 1 and com- 
prises two separate sections: one that performs fade identification and the other that per- 
forms fade prediction. In the identification section, D1 and D2 perform as discriminators 
carrying out the operations A obs {t ) - A {x A (f,))J and | A obs (t ) - 2 A (r, )jj, respec- 

tively, multipliers Ml and M2 form the operations needed to complete the emulation of the 

relations given in Eq.(23), the inverters I perform multiplicative inversion, and 5 denotes 
the emulation of the operation yielding A m exp {^i nA x A (t )J . Also, denotes the multi- 
plication by the discrete time transition factor exp (— y s T- j , &. is the square of this factor, 

and 7/ is a time delay of 7/ seconds. In the prediction section, once a time interval T pre d 
is selected over which a fade prediction is to be made, D (T pred ) evaluates the operation 
1 — exp ( - 2 y s T pred ), y pred then performs the operation as specified by Eq.(34), and fi- 
nally multiplier M3 forms the second member of the right side of Eq.(33). 

Such an algorithm can be implemented either in software or hardware and is easily 
included into the operation of any satellite earth station. Although the structure of the 
algorithm remains the same in every application, several of the numerical coeficients which 
it employs reflects the specific geometry and operating frequency of a satellite link as well 
as the location of the earth terminal. For example, y s is a function of the propagation path 

length through the potential rain region as well as the elevation angle, as shown by Eq.(4). 
Not only is the link operating frequency reflected in the selection of the parameters A m and 
<T /n/1 that appear in each 5 module of Fig. 1, but also the specific location of the earth sta- 
tion due to the extremely detailed rain statistics data base of the ACTS rain attenuation 
model (a geographical resolution of ~ 0.5° in latitude and longitude in the continental US, 
essentially equivalent to about 6000 rain "zones") from which they are calculated. The 
software implementation of the ACTS Rain Attenuation Prediction Model for this and other 
general satellite system link design is available from COSMIC, NASA Software for Indus- 
try, The University of Georgia, Athens, GA 30602, U.S.A. 
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APPENDIX - A Derivation of Equation 13 


Here, it is desired to augment the dynamic description of the random process x A ( t ) 

that is afforded by the Kolmogorov equation, Eq. (10), with observed a priori measure- 
ments A 0 bs(ti) of the link attenuation which, through Eq.(8), is a convolution of x A (f, ) as 
well as the measurement noise n (r). It is thus expedient to consider a quantity that incor- 
porates both aspects of a priori observation and dynamical transition, i.e., the conditional 
transition probability density p (t +At ),x A (t +At),A o6j (r t - +At )l A ofcj (r, )J which 

incorporates two successive observations at time t i and at time t i +At with At = t i+l - t i 

used to explicitly denote the prevailing time increment, and with the proviso that 
t i- t<t i+ 1 . This probability density can be decomposed into two equivalent cascaded 
forms 

P ( A obs (fi + At )l A obs (ti ) , x A (t + At ) , x A (t + At )) x 

x p ( {t + At ), x A (t + At )l A obs (ti)) 
and 

P f ■ x a(* +a ' )>* A (t +Ar )l A ob M,A obs (ti +At ))p ( A +At )l A 

Equating these two representations and rearranging factors gives 

p\x A {t +At ),x A (t +At )l A obs (ti), A obs (ti +At )) = F(At ) x 
x P ( A obs (ti + At )\ A obs (ti),x A (t + At ),x A ( t + At )J x 

x p ( jr A (t + At ) , x A (t + At )l A obs (ti )) (A 1 ) 

where F (At ) = J A obs (t t + At )l A ofci (r,)j is a factor that is independent of the process 
x A . The first probability density function on the right side of Eq.(Al) is connected with the 
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uncertainty associated with the measurement process. Noting that the measurement noise at 
+ At only affects the random quantity A obs {t i + At) and since it is only the sampled 
value x A (t + At) that is connected with this sampled Gaussian noise process as defined by 
Eq.(7), one has the following development: 


P A obs (ti + At )l A obs (ti), x A (t + At ) , x A (t + At )j = p ( A obs {ti + At )\x A (t + At) 


= N exp 


n(ti + At) 


2<r„ 


At 


’[♦( 


= N exp\<p[x A (t + At),t+ At 


) 4 <] 


N 1 + <p{x A (t+ At),t+ At ^At +... 


(A2) 


where N is a normalization constant and use was made of Eqs.(12) and (14) and the fact 
that At is sufficiently small to warrant the series expansion in the last line. The second 
probability density in Eq.(Al) is an extrapolated probability density function of the process 
x A based on the a priori observation A obs (t i ). Expanding this quantity into the first two 
terms of its series representation yields 

P^x A (t + At), x A (t + At )\ A obs (tj)J = p(x A (t),x A (t )l A obs (t i )) + 

dp(x A (t),x A (t)\A obs (t i j) 

+ 1 ’-At (A3) 


After the initial a priori input A obs (ti) at time //, no other such information is used during 
the time interval At thus allowing the time derivative in to be described by the Kolmogorov 
equation, Eq.(10). Making this substitution and inserting this result and that of Eq.(A2) 
into Eq.(Al) and retaining terms up to first order in At yields 


P ( x A (r +At ),x A (t +At )l A obs (ti ) , A o6j (t, +At )j = C 


p(x A (t),x A (t)\A obs (t i )) 


+ 
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+ D x\p ( *A (* )> *A (' )! A ofcs(*i)j| At + 

+ <p (x A ( t +At ),(t +At )) p ( X A (t ) , x A ( t )l A obs {ti )) At 


(A4) 


where C = F(At)N is a constant (at least for a fixed At ) which must now be evaluated. To 
this end, one integrates both sides of this equation over all values of x A (r ) and x A ( t ) and, 
noting the analogous relationships that follow Eq.(15), obtains 


C = 


-.-l 


x A ( t +At ),(t +At )) p ( x A (t) , x A (t )l A )) dx A dx A 


» 1- At 


U" 

J j <t> (t +At }>(' +At)}p(x A (t),x A (t )l A (q )) dx A dx A (A5) 


where, once again, the fact was used whereby At is small. Substituting this result back 
into Eq.(A4) and retaining terms up to first order in At yields 


p ( x A (t +At ) , x A (t +At )l A obs (t t ),A ^(ti +At )) = 


= p ( * A (t ) , x A (t )\A obs {ti)) +D Xa p(x A {t),x A {t)\A oi(S (r ( ))j At + 


+ 


<P^ A (t+At),(t +At 



X 


P ( ) , x A \t )l A obs (ti )) dx A 'dx A ' 


p(x A {t),x A {t )l A obs (ti))At 


Finally, moving the first member on the right side to the left side and dividing through by 
At forms a different quotient which reduces to Eq.(13) after taking the limit At — >0. 
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